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ABSTRACT 



We propose a physically motivated and self-consistent prescription for the modeling 
of transient neutron star (NS) low-mass X-ray binary (LMXB) properties, such as duty 
cycle (DC), outburst duration and recurrence time. We apply th is prescription to th e 



population synthesis (PS) models of field LMXBs presented by iFragos et al.1 (|2008l ). 
and compare the transient LMXB population to the Chandra X-ray survey of the two 
elliptical galaxies NGC 3379 a nd NGC 4278, which revealed several transient sources 



(jBrassington et al.l 120081 . |2009j). We are able to exclude models with a constant DC 
for all transient systems, while models with a variable DC based on the properties of 
each system are consistent with the observed transient populations. We predict that 
the majority of the observed transient sources in these two galaxies are LMXBs with 
red giant donors. Our comparison suggests that LMXBs formed through evolution of 
primordial field binaries are dominant in globular cluster (GC) poor elliptical galaxies, 
while they still have a significant contribution in GC rich ones. 
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1. INTRODUCTION 



Recent Chandra observations (|Kim et al.l 12006 : iBrassington et al.l 120081 . 12009 ) have yielded 
the first low-luminosity X-ray luminosity functions (XLF) of LMXBs for two typical old elliptical 
galaxies, NGC 3379 and NGC 4278. The detection limit in these observations was ~ 3 x 10 36 erg s _1 , 
an order of magnitude lower than in mos t previous similar s urvey s , whil e the completeness limit 



was ~ 10 37 ergs 1 . In a followup study, IBrassington et ajj (|2008l . [2009) used multiple Chandra 



pointings to identify potential transient sources. In that work, a transient candidate was defined as 
a source with X-ray luminosities varying by at least a factor of 10, above the completeness limit, 
between different pointings. Similarly a potential transient candidate was defined as a source with 
X-ray luminosities varying by at least a factor of 5. Using these definitions, 5 transient candidates 
and 3 additional potential transient candidates were identified in NGC 3379, while for NGC 4278 
the corresponding numbers were 3 and 3 respectively. 



Piro Bildstenl (|2002l ) conjectured that LMXBs in the field of ellipticals must be dominated by 



NS accreting from red giants. Using analytical approximating formulae for the evolution of binary 
systems they estimated that such binaries are expected to be transient, with maximum recurrence 
times of 100 — 10000 yr, for at least 75% of their lifetimes. However, they limited their study only 
to LMXBs with red giant donors, and were unable to calculate the duration of the outburst phase 
or to predict how many transients would be identified in repeated Chandra visits. 



Fragos et al.l (|2008l ) presented results from extensive LMXB population synthesis (PS) sim- 



ulations for NGC 3379 and NGC 4278. They considered models for the formation and evolution 
of LMXBs from primordial and isolated field binaries with different common envelope efficiencies, 
stellar wind prescriptions, magnetic braking laws, and initial mass functions. They identified mod- 
els that produce XLFs consistent with observations both in shape and normalization, suggesting 
that a primordial galactic field LMXB population can have a significant contribution to the total 
population of an elliptical galaxy. 



In this letter we expand on the work bv lFragos et al.l (|2008l ). and propose a physically motivated 
and self-consistent prescription for the modeling of transient LMXB properties. Furthermore, we 
carry out Monte Carlo s imula tions to compare our findings to the recent observational work by 
Brassington et al.l (|2008l . |2009| ). In § 2 we describe briefly the PS models we use for the formation 
and evolution of field LMXBs . In § 3 we expand on the detailed modeling of transient LMXBs 
and our comparison of PS models with observations. We show our results in § 4, and discuss the 
implication of our findings for transient LMXBs in § 5. 
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LMXB POPULATION SYNTHESIS MODELS 



The models presented in the study by iFragos et al.l (|2008l ) were focused on LMXBs formed 
in the galactic field as products of the evolution of iso l ated p rimordial binaries. The simulations 
were performed using StarTrack (jBelczvnski et al.ll2002l . 120081 ). an advanced PS code that has been 
tested against and calibrated using detailed MT calculations with stellar structure and evolution 
codes, and observations of binary populations. In the development of these models, all current 
knowledge about t he stellar population in these galaxies were incorporated (see Tables 1 and 2 in 
Fragos et al.ll2008l ). The observationally determined parameters, such as their age (9.5 — 10.5 Gyr) 
and metallicity ([Fe/H] = 0.14 - 0.16), or their total stellar mass (8.6 x 10 10 - 9.4 x 10 10 M ), are 
similar for NGC3379 and NGC4278. However, there are physical processes involved in the formation 
and evolution of a binary system, such as stellar winds, initial mass function, common envelope 
efficiency, magnetic braking and transient behavior, which are not fully understood quantitatively. 
Therefore various prescriptions were used to model them. A total of 336 models were studied. See 
Tables 3 and 4 in lFragos et al.l (|2008l ) for all the model parameters and naming conventions. 



Fragos et al.l (120081 ) found that a small subset of their models produce XLFs in very good 



agreement with the observations, based on both the XLF shape and absolute normalization. They 
concluded that formation of LMXBs in the galactic field via evolution of primordial binaries can 
have a significant contributio n to the total populat ion of an elliptical galaxy, especially the GC 
poor ones such as NGC3379 (jFabbiano et al.l 120071 ). In our analysis here, we adopt their best- 
fitting model (hereafter model 14 /T ), which produces the XLF in best agreement with observations. 
Model 1A IT assumes a moderate commo n envelope efficienc y acE = 0.5, a standard stellar wind 
prescription pres c riptio n (as described in iHurley et al.ll2002l ). the magnetic braking law derived by 
Ivanova &: Taaml (120031 ). and a Salpeter initial mass function. 



The LMXB population predicted by model 14^ T has a significant contribution from transient 
systems, which with reasonable outburst DC can even dominate the XLF. A s a consequence the 
XLF shape is rather sensitive to the treatment of these transient systems. IFragos et al.l ((2008) 
showed that keeping the same PS parameters and changing only the modeling of transient sources 
leads to different XLFs. They tried different methods of modeling the outburst characteristics of 
transient LMXBs and got the best agreement with observations when they consider a variable DC 
for NS LMXB based on the binary properties, hereafter model 14A IT . However, a constant low DC 
(~ 1%), hereafter model 1AB IT , was also statistically consistent with the observed XLF. In both 
cases, the outburst luminosity was calculated as described in section 3.1. 
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3. MODELING THE TRANSIENT BEHAVIOR OF LMXBS 

3.1. Duty Cycle of Transient X-ray Binaries 

In our models we keep track of all binary properties, including the MT rates, as a function 
of time for all accreting compa ct objects. Here we focus only on NS LMXBs, which dominate the 
total population (see Fig. 2 in iFragos et al.ll2008l ). We use the MT rate M d to identify transient 
sources in our simulation results. In the context of the thermal disk instability model, mass trans- 
ferring b inaries with MT rate lower than the critical rate M rr jt for the thermal d isk instability to 



develop (Ivan Paradiji 



1996 



King et al 



1996 



Dubus et al. 



1999 



Menou et al.ll2002l ). are considered 



transient sources, meaning that they spend most of their life in quiescence (T qu i eS cence ) , when they 
are too faint to be detectable, and they occasionally go into an outburst. The fraction of the time 
that they are in outburst (T out b U rst) defines their DC: 



DC 



1 outburst 



/(Tc 



outburst 



quiescence j 



(1) 



The details of the thermal disk instability model are not well understood therefore the values 
of the outburst luminosity and the DC cannot be calculated directly from first principles. A simple 
but physically motivated treatment is to assume that in the quiescent state the NS does not accrete 
any (or accretes an insignificant amount of) mass, and matter from the donor is accumulated in 
the disk. In the outburst state all this matter is accreted onto the NS emptying the disk. Taking 
into account also that t he X-ray luminosity probably cannot exceed -^Edd by more than a factor of 



2 (cf. iTaam et al.lll997l ). we define the outburst luminosity as: 



T -for GM a M d 1 \ 

L x = %oie x mm I 2 x L Edd , — — x — I , (2) 

where M a and R & are the mass and the radius of the accretor, dotM& is the mass-loss rate of the 
donor and r/bol is correction factor that converts the bolometric luminosity to the observed Chandra 
band. 



In equation [21 DC is unknown. iDobrotka et al.l (|2006l ) studied accretion disk models for dwarf 



novae which are thought to experience the same thermal disk instability. They found a correlation 
between the DC of the system and the rate at which the donor star is losing mass M d . The 
exact relation between these quantities depends on the disk's viscosity parameters, but the general 
behavior can be approximated by: 

DC^(M d /M cr S . (3) 



Plugging equation ([3]) into equation (J2j) we eliminate the DC dependence and get an expression 
that depends only on quantities which are directly calculated in our PS modeling. 
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3.2. Recurrence Time and Duration of X-ray Outbursts 



The derivation of the critical MT rate M cr ; t for the thermal disk instability, comes from re- 
quiring that some part of the accretion disc has a temperature below the hydrogen ionization 
temperature. In the thin disk approximation the effective temperature decreases with increasing 
disk radius, so it is sufficient to require that the effective temperature at the outer edge of the disk is 
below the hydrogen ionization temperature. In X-ray binaries wit h NS accretors, the surface tern 



perat ure of the disk is regulated by irradiation from the accretor (jvan Paradiis 
200lh . 



1996 



Dubus et al 



Following the analysis presented bv lPiro Bildstenl (|2002l ). we assume that the accretion disk 
around the NS extends up to a radius of 70% of the Roche lobe radius of the NS (-Rrl.i), and that 
a transient NS LMXB is in the quiescent state until the surface density at the edge of the disk 
reaches the critical surface density 



644(Mi/M ( 



-0.37 



(R/Rq) 1 - 11 gcm~ 2 



(4) 

In equation [J] we assume a pure hydrogen disk 



that will lead to instability (jDubus et al.l 120011 1 
with viscosity parameter a = 0.1, and no irradiation of the disk while the binary is in the quiescent 
state. For systems with helium accretion disks , such as LMXB w ith white dwarf donors, we 
modify equation accordingly (see equation 2 in lLasota et al.ll2008l ). Hence, the maximum mass 
accumulated in the accretion disk, before the X-ray binary goes into outburst, is: 



M, 



disk, max 



0.7R m 



Rns 



T, msx {R)2^RdR 



(5) 



From the equation above we can set an upper limit on the duration of the quiescent phase and the 
outburst. Using the mass-loss rate of the secondary star, as given by our PS models, and equations 
[T] and [3l the time duration of the quiescent and outburst phase becomes: 



-^disk,: 



- quiescence 



M d 



f outburst 



Ml 



M c 2 rit 



Ml 



-M, 



disk, max- 



(6) 



3.3. Detection probability of transient LMXBs 

The probability that a transient LMXB is detected as a bright X-ray source is equal to its DC. 
However, the probability that a LMXB is identified as a transient source in a survey with multiple 
observations depends on the duration and the recurrence time of the outbursts of the source, as 
well as the time spacing between the observations. In our analysis we assume that a transient 
LMXB can be in two states, outburst or quiescent, and that its luminosity in outburst is given by 
equation [21 while in quiescence it is undetectable (i.e. we do not take into account rise and decay 
characteristics of the outburst). Hence, in order for a source to be identified as a transient, it has 
to be found at a different state, in at least one out of the multiple observations (see Figured]). 
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Fig. 1 . — Schematic description of the procedure followed to calculate the probability that a LMXB 
is identified as a transient source in a survey. 
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Given a model population of LMXBs, we calculate for each transient LMXB the probability 
that it would be identified as a transient s ource in an obse rvational survey like the ones for the 



elliptical galaxies NGC 3379 and NGC4278 (|Kim et al.ll2006l ). In order to estimate this probability. 



we first calculate the outburst and recurrence time for the specific binary as described in section 
3.2, and then create an artificial lightcurve. This lightcurve should include at least one full cycle 
of the transient and be at least twice as long as the time interval between the first and the last 
observation. Next, we choose randomly a time in this artificial light curve to position the first 
pointing of the survey with the other pointings following according to the time intervals between 
the pointing of the survey we are simulating. For this random positioning of the first pointing, 
we attain the state of the transient LMXB (outburst or quiescent) at each of the pointings, and 
assess whether this system would be identified as a transient source (see Figure [1]). We repeat 
the procedure 10 5 times, choosing randomly a different starting time in the artificial lightcurve, 
which enable us to estimate the probability that a specific system would be identified as a transient 
source. 

Based on this probability for each transient LMXB of the model population, we perform 
another Monte Carlo simulation, this time to estimate the total number of transients that would 
be identified in the modeled survey. Repeating this procedure 10 5 times, we are able to derive a 
probability density function (PDF) of the number of identified transient sources, for a combination 
of PS and X-ray survey parameters. 



4. RESULTS 



As mentioned earlier, we adopt the PS model 14 /T , which was chosen b ased on comparison 
of the shape of the modeled XLF to the observed one (IFragos et al.l l2008l ). We examine two 
prescriptions for the transient characteristics of LMXBs: a variable DC based on the properties 
of each binary as described in section 3.1 (model 14^4 /T ), and a low consta nt DC of 1% (mode l 
UA IT ). For more details on the model parameters see sections 2.3 and 3.1 in lFragos et al.l (|2008l ). 
We normalize the total number of LMXBs in our models in two different ways: based on the total 
stellar mass of the galaxy, or based on the number of observed LMXBs above the completeness 
limit of the observations. 

In order to simulate the observational identification of transient sources as realistically as 
possible, for each of the two elliptical galaxies we use the actual time intervals between the successive 
pointings of the survey. Furthe rmore, we try to m i mic the defi nitions of transient candidate and 
potential transient candidate of iBrassington et al.l (|2008l . |2009| ). In our simulated observations, a 
transient candidate is a system that is identified as transient and has outburst luminosity above 
10 38 ergs -1 , while a potential transient candidate has outburst luminosity between 5 x 10 ergs -1 
and 10 38 ergs -1 . 



Figure [2] shows the PDF for the number of detected sources in an X-ray survey such as that 
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Fig. 2. — PDF o f the numbe r of detected transient sources in an X-ray survey such as that by 
Brassington et alj ( 20081 . 12009I ). for model 14A IT in panels (a) and (b) and model 14B IT in panels 



(c) and (d). Panels (a) and (c) correspond to the elliptical galaxy NGC3379, while panels (b) 
and (d) refer to NGC4278. Solid lines are used when the normalization of the total number of 
LMXBs in our models is based on the total stellar mass of the galaxy, whereas dotted lines signify 
a normalization based on the observed number of LMXBs above the detection limit. The vertical 



dashed lines refer to the actual observations by iBrassington et alJ (|2008l . I2009I ) 
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by iBrassington et al.l (|2008l . |2009| ). When we normalize the total number of LMXBs to the stellar 
mass of the the galaxy, model 14A IT (variable DC) gives PDFs consistent with observations for 
both transient candidates and potential transient candidates, and for both galaxies. On the other 
hand, model 14B IT (constant DC ~ 1%) is consistent with observations only when we look at the 
bright transient sources (i.e., transient candidates). A constant DC significantly overproduces the 
number of potential transient candidates. 

Normalizing the total number of LMXBs of our models to the observed number of LMXBs 
above the completeness limit of the observation (~ 10 37 ergs -1 ), for each of the two galaxies 
separately, alters our findings. In the case of NGC3379, and for model 14A IT , the PDFs of 
both transient candidates and potential transient candidates are now in excellent agreement with 
the observation s. This suggest t hat our models slightly overproduce LMXBs. This was already 



mentioned in iFragos et al.1 ([2008), where the authors explained that the PS model normalization 
can be very sensitive to certain PS parameters, such as the distribution of the mass ratios between 
the initial masses of the two binary components. Without a proper multi-dimensional coverage of 
the parameter space of at least thousands of models and given the uncertainties in the total galaxy 
mass, we cannot use the total number of LMXBs as a formal constraint to a factor better than a 
few. For the case of NGC 4278 on the other hand, normalizing our model population to the number 
of observed LMXBs results to a very high number of expected identified transient source, which is 
inconsistent with observations. However, this does not come as a surprise, since the elliptical galaxy 
NGC 4278 has a large number of GC (~ 5 times more than NGC 3379) and a significant fraction 
of the observed LMXB population resides in GC. It has already been suggested that lum inous GC 
LMXBs are predominantly persistent source (jBildsten fc Deloydl2004l : llvanova et al.ll2008l ) in which 
case the similar number of identified transients in the two elliptical galaxies, despite the difference 
by a factor of ~ 2 in the total number of observed LMXBs, is expected. 



Fragos et al.l (|2008l ) analyzed their model LMXB populations and found that the sub-populations 



that mainly contribute to the model XLFs are transient and persistent NS-LMXBs with red giant 
donors. Following a similar analysis, in Figure[3]we show again the PDF of the number of potential 
transient candidates for models 14A IT and 14B IT of the elliptical galaxy NGC 3379, but this time 
we split the modeled population into sub-populations based on the type of the donor star. We 
find that in model 14A IT , the majority of i dentified transien t sources are LMXBs with red giant 



dono rs, as already argued in the literature (jKing et al.1 Il996l : iPiro &: Bildstenl [2002J; IFragos et al 



20081 ). In contrast, model 14B IT predicts a much higher than the observed number of identified 
transient LMXBs and the majority of them have main sequence donors. This results from assigning 
the same low DC to all systems, whereas transient LMXBs with main sequence donors have usually 
MT rates close to M cr i t . These systems would otherwise have DC ~ 1 and low outburst luminosity, 
and thus would be undetectable at X-ray luminosities above 5 x 10 37 ergs -1 . 
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DISCUSSION 



The goal of t his work was to po s e further co nstraints on our models based on the recent obser- 
vational work by iBrassington et al.l (|2008l . 120091 ). and simultaneously gain a better understanding 
on the nature of transient LMXBs. We foun d that the populations of transient LMXBs produced 
by the PS model 1AA IT ( Fragos et al.ll2008l ) is consistent with the observed population. In this 
model a variable DC was assigned for each system based on its properties (see section 3.2). In 
contrast, we show that the widespread in the literature assumption of a constant lo w [~ 1%) D C 



for a ll transients, although it can result in a XLF consistent with the observed one (IFragos et al 



20081 ). is not consistent with the number of observationally identified transient LMXBs. 



While the two elliptical galaxies modeled here have very similar properties in terms of stellar 
mass, age, star formation history and metallicity, the major difference between the two is their 
specific GC frequ ency. NGC 4278 has about 5 times more GC per unit mass compared to NGC 3379. 
Kim et al.l (|2006l ) using similar Chandra observing times for both galaxies, found that NGC 4278 



has ~ 3 times more LMXBs with luminosity above 10 37 ergs _1 than NGC 3379. The excess of 
X-ray sources in NGC 4278 can be attributed to the higher LMXB formation efficiency in GC, 
suggesting that for the GC-rich elliptical galaxies the formation of LMXBs in GC through dynamical 
interactions plays an important role. At the same time, the number of identified transient LMXBs 
is roughly the same in both elliptical galaxies. Despite the small statistical sample of only two 
galaxies, it seems that the number of transient sources is proportional to the stellar mass of the 
galaxy rather than the number of its GC. 

The same argument is also supported by our PS models which only take into account LMXB 
formation through binary evolution of field primordial binaries. When we normalize the modeled 
LMXB population to either the stellar mass of the galaxy or to the number of observed LMXBs 
in the GC-poor galaxy NGC 3379, the predicted numb er of observable transient LMXBs is consis- 
tent with observations (jBrassington et al.l 120081 . |2009| ). Hence, we conclude that LMXBs formed 
through the evolution of primordial field binaries dominate the observed LMXB populations in 
GC-poor elliptical galaxies, while in GC-rich ellipticals the formation of LMXB through dynamical 
interactions in GC becomes more important. However, even in the latter case, the primordial field 
population has a significant contribution. 



This work was supported by Chandra GO grant G06-7079A and subcontract G06-7079B. 
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Fig. 3. — PDF of the number of potential transient candidates for the model 1AA IT (variable 
DC) and model 1AB IT (constant DC). We show the contribution of different sub-populations by 
separating the LMXBs into groups of systems with different donor. The vertical dashed lines refer 
to the observed number. 



